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We predict that guiding center (GC) diffusion yields a linear and non-saturating (transverse) 
magnetoresistance in 3D metals. Our theory is semi-classical and applies in the regime where the 
transport time is much greater than the cyclotron period, and for weak disorder potentials which 
are slowly varying on a length scale much greater than the cyclotron radius. Under these conditions, 
orbits with small momenta along magnetic field B are squeezed and dominate the transverse conduc¬ 
tivity. When disorder potentials are stronger than the Debye frequency, linear magnetoresistance is 
predicted to survive up to room temperature and beyond. We argue that magnetoresistance from 
GC diffusion explains the recently observed giant linear magnetoresistance in 3D Dirac materials. 


Magnetoresistance provides a powerful means with 
which to probe the scattering history of particles in 
a magnetic field. Departure from the conventional 
paradigm - quadratic magneto-resistance at low fields, 
saturating at high fields [ 1 ] - signals anomalous particle 
scattering behavior. One particularly appealing regime 
is non-saturating and linear magnetoresistance (LMR), 
which has a long standing history mg given its poten¬ 
tially disruptive technological impact [ 6 ] . 

Very few theories predict LMR in a closed single com¬ 
ponent Fermi surface. A well known example is Ref. [7], 
which showed that Dirac metals in the extreme quantum 
limit (when only the n = 0 Landau level is occupied) 
exhibit LMR in the presence of screened Coulomb im¬ 
purities. Another mechanism yielding quasi-linear MR 
arises from inhomogeneity [HHin]- However, significant 
LMR in these requires strong inhomogeneity [10]. Con¬ 
temporary proposals that extend the above treatments, 
have also found LMR under similar requirements [TTJ |T2] . 

Recently, giant LMR that lie outside the above two 
paradigms [see (i) and (ii) below] was reported in the 
newly discovered class of three-dimensional Dirac mate¬ 
rials (3DDM) [TBUTT] . LMR in 3DDM exhibit puzzling 
features including (i) its occurrence when multiple Lan¬ 
dau levels are occupied far from the extreme quantum 
limit, and (ii) arising in weakly disordered, high mobil¬ 
ity samples. Further, LMR manifests consistently over a 
variety of 3DDM experiments, including in TiBiSSe m, 
Cd 3 As 2 [Hills], Na 3 Bi[T 6 ], and TaAs m, where chemi¬ 
cal potential /i typically lies 0.1 eV above the Dirac point, 
hinting at a single underlying explanation. 

Here we propose a semi-classical mechanism for LMR 
in metals, wherein charge transport is dominated by 
guiding center (GC) motion. Importantly, this mecha¬ 
nism naturally gives giant LMR under (i) and (ii) above, 
explaining the puzzling behavior [IMg. The main re¬ 
quirement is that the disorder potential is smoothly vary¬ 
ing on a scale, which is large as compared to the cy¬ 
clotron radius, re- The main features of GC magnetore- 
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FIG. 1: a) Magnetoresistance can be dominated by guiding 
center (GC) motion when disorder correlation lengths ^ ^ rc. 
Here rc is the cyclotron radius. This regime is characterized 
by slow GC motion, Vgc = VrU(r) x z/B, accompanied by 
fast cyclotron orbits, Veyd- GC diffusion in this environment 
givers rise to LMR. b) Electrons perform closed orbits of the 
Fermi surface, with GC motion classified into two types, kz < 
k* (red), and kz > k* (blue); critical kl (green), c) For 
kz <kl, electrons are squeezed in z yielding mean free paths 
I ^ and in-plane Dxx ~ 'Cgc^ (x 1/B. d) In contrast, kz > 
k* electrons exhibit unconstrained motion yielding in-plane 
Dxx ~ vlc^/vz oc 1/B^ (see text). 


sistance are exposed by writing the transverse resistivity 
as 


_ ^xx _ Q ^ _ tan^^ , . 

^xy l + [tan6>//]^ 

where cr^x and cr^y are the transverse {x-y plane) conduc¬ 
tivity and Hall conductivities respectively, and tcHiOH = 
^xy!^xx is the Hall angle. Using the familiar — nejB 
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with n the density, and e the carrier charge, we have 


As we argue below, in the regime of cjc^tr ^ 1 and ^ ^ 
Tc, GC diffusion gives a Hall angle, and therefore that 
is independent of magnetic field magnitude, leading to 
LMR in Eq. Q. Here uJc is the cyclotron frequency, Ttr 
is the transport time, and the magnetic field B = Hz. 

Guiding center magnetoresistance can be understood 
as follows. In semi-classically large B fields (cJcRr ^ 1)? 
electrons exhibit in-plane trajectories (t) characterized 
by slow guiding center motion R(t) accompanied by fast 
cyclotron orbits rcyci(t). The latter, characterized by Tc, 
depends on intrinsic material properties and B; whereas 
the former depends on the potential profile sampled by 
the electron over one cycle which can yield unusual tra¬ 
jectories m- A unique situation arises for slowly varying 
disorder potentials H(r). In this regime ^ ^ Tc (see Fig. 
la), electron trajectories are dominated by guiding cen¬ 
ter motion which follows the local disorder landscape at 
H, with velocity Vgc = [VrE(R)] x zfB. 

Guiding center diffusion is characterized by diffusion 
constant ^ The central question is: what is r? 

First, it is important to note this picture is not valid in 
strictly 2 D, because GGs form closed orbits along equipo- 
tential lines. Hence it is crucial to include motion in the 2 : 
direction which restores diffusive motion. There are two 
classes of electron motion depending on their kz value 
with respect to (Fig. lb). For kz > k^^ electrons 
possess kinetic energy in the z direction exceeding the 
typical potential fluctuation. As a result, the electron 
moves freely across many potential fluctuations shown in 
Fig. Id. Within time r> ^/vz^ the GC senses a differ¬ 
ent local electric field and changes direction. As a result, 
> kl) ^ ^lc^> oc 1/H^. Using a^y = ne/B and 
Eq. we recover the standard saturated magnetoresis¬ 
tance. 

On the other hand, electrons with kz < k^ are typ¬ 
ically squeezed by a local potential barrier. As shown 
in Fig. Ic, they must travel sideways by a distance ^ 
to get around the barrier. In this case, r< « 
and ^ ^ IjB; LMR follows immediately from 

Eq. § and Q. Paradoxically for large H, it is the elec¬ 
trons squeezed in the z direction which dominate trans¬ 
port in the x-y plane, leading to LMR. The importance of 
squeezed electrons was pointed out by Ref. muo] in the 
context of magneto-transport of Boltzmann gases. We 
adapt their reasoning to the case of a degenerate Fermi 
sea and present a more quantitative treatment below. 

We emphasize that our mechanism for LMR does not 
have anything to do with the Dirac spectrum per se. Nev¬ 
ertheless, 3D Dirac semi-metals provide an ideal venue 
that satisfy the conditions required for LMR. First, the 
high mobility of 3DDM (77 > few x lO^cm^/Vs) allows 
the regime where electrons undergoes many cyclotron or¬ 


bits before scattering, cJcTtr = By ^ 1, to be achieved 
at relatively low magnetic fields. Second, the relatively 
small chemical potential y ^ 100 meV but large Fermi 
velocity vp 1 — 10 x 10 ^ cm /sng, give small cyclotron 
radius Vc ^ ji/evpB = 10 — lOOnm even at IT. Third, 
large dielectric constants of ^ 40 [2TJ |22] effectively 
screen Coulomb impurities to yield weak and slowly vary¬ 
ing disorder potentials, with large correlation lengths, 
^ ^ 20 — 60 nm [22]. As a result, ^ ^ Tc at relatively 
low H ^ 1 T, allowing GC motion to dominate the x-y 
plane (transverse) magnetoresistance. 

We note parenthetically that our regime of interest is 
distinct from multi-component systems e.g., charge com¬ 
pensated systems which can be ultra-sensitive to mag¬ 
netic field [23[ |24|. Instead we are interested in LMR 
in 3DDM [T3IIT7] which were observed in the metallic 
regime with carriers in a single band [14]. 

We begin by considering the diffusive motion of 
charged particles in a magnetic field, B = Hz, and a 
slowly varying and weak disorder potential V (r) . While 
formally interested in 3DDM, our analysis below is gen¬ 
eral; we will only specify 3DDM as needed to com¬ 
pare to recent experiments. Disorder is characterized 
by (U(r)U(r')) = VqB{\y — r'|/^), where {O) denotes 
disorder averaging, and ^ the correlation length. T is 
a dimensionless function that vanishes for |r — r'| ^ 
Lastly, we will be interested in weak disorder strength 
eVb < /i seen in 3DDM experiments [IMZl, and recent 
estimates [22] . 

Equations of motion - The motion of particles on the 
Fermi surface with chemical potential y can be described 
by the semi-classical equations of motion 

mv_L = —eVrU(r) + ev_L x B, (3a) 

mvz = -edzV{T), (3b) 

where m is the cyclotron mass, v±{kz, y) = {vx^Vy), and 
Vz{kz^y) are velocities transverse to the magnetic field 
and along the magnetic field respectively. We note that 
throughout our analysis below, these quantities depend 
on momentum along the field, kz, and y. For e.g., ve¬ 
locity is captured via group velocity Vk = h~^ds\i/d]<. so 
that the x-y plane speed for Dirac particles is | | = '^0 = 

VF\/l — k\/k\, and m = yjv^ [25], where vp is Fermi 
velocity, and kp the Fermi wave vector; 5k is the particle 
dispersion. For brevity, we will drop explicit mention of 
kz dependence, bringing it up when necessary. We do 
not expect Berry phase related terms [26] to contribute 
to the transverse magnetoresistance behavior that we are 
interested in here. 

The trajectories of charged particles, r(t), in crossed 
B and V (r) can be complex, since they involve transport 
processes spanning multiple time scales (e.g., cyclotron 
period, guiding center scattering time, and transport 
time). However, in semi-classically strong fields {oOcXtr ^ 
1), and for a slowly varying potential so that correlation 
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length is larger than cyclotron radius rc), its motion 
is conveniently captured via r(t) = R(t) + rcyci(t). Here 
R(t) is the slow moving 3D guiding center coordinate, 
whereas rcyci(t) describes fast cyclotron motion lying in 
the x-y plane. 

This reasoning yields the following ansatz for velocity 
in the = {rx^Vy) plane as m 


We adopt a mean-field approach in estimating r. Since 
T rapidly decays for AR > ^, r is most sensitive to the 
way the GC moves in AR{t) < ^. As a result, we write 
dAR = with speed 'r’av = [('^gc)| + aver¬ 

aged over a single domain; here (O)^ denotes averaging 
across a single domain. Changing variables t AR, 
yields 


v±{t) = -h i)gc(^), 


lygc 


iE{r^) 
B ’ 


( 4 ) 


where ujc = eBjira^ and used complex notation Ox + 
iOy = O for vectors in the x-y plane. The latter part 
of Eq. Q was obtained by substituting the ansatz into 
Eq. (3r) and setting mhgc = 0 for slowly varying V{r). 
Eq. (4) is valid for |mhgc| <C \eE{r^)\. Estimating E « 
Vb/C^ we obtain the condition 


eW 




^0 


^where Vc = ^. (5) 


UJc 


Since we are interested in weak disorder eVb < the 
above condition is satisfied within our regime of validity, 
^ > rc. 

Motion in 2 : can be understood in the following way. 
Eirst, we note that for ^ ^ Vc and uJcBy ^ 1, the po¬ 
tential the electron feels is determined by (r(t))i cycle = 
R(t). Next, for Vgc ^ Vz electrons, the GC moves slowly 
in the x-y plane as compared with z. As a result, inte¬ 
grating Eq. yields energy conservation 

^ _e{y[R(t)]-l/[R(0)]}, ( 6 ) 


where we have set Vr^ V (r) •5tr_L = 0. This is valid when 
|Vr^E(r) • dtrj_\ < \vzdzV{r)\. Estimating \dtr±\ « v^c 
and using disorder that is isotropic, yields the original 
condition Vgc ^ Vz. 

Guiding Center Transport - The separation of time 
scales between slow GC motion, and fast cyclotron mo¬ 
tion enables us to write the velocity correlator as 

{v±{tWm « (7) 


where we have used a relaxation-time approximation in 
the last term to capture the Drude contribution to mag¬ 
netotransport [18]. 

Replacing r(t) with its average over one cycle as 
above, and using Eq. Q, we find GC diffusion, = 
( 1 / 2 ) /o“ {vsc{t)vgc{0))dt, as 

pOO 

/ {Em)]E[B.m)dt/{2B^)=ElT/{2B% 

JQ 

dtE{AR/i), (8) 

where AR = |R(t) — R(0)|, Eq is the characteristic elec¬ 
tric field strength of the disorder potential, and r is the 
scattering time that is sensitive to the GC trajectory. 
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( 9 ) 


where M is a number of order unity. Using gaussian cor¬ 
relations, (U(x)U(O)) = VqE{x) = VQe ~^^obtain 


A = v^/ 2 , and Eq(^ = 6Vq. 

Two distinct classes of GC trajectories can be dis¬ 
cerned: (a) squeezed z-motion (Eig. Ic), and (b) unre¬ 
stricted z-motion (Eig. Id). Squeezing in class (a) arises 
from energy conservation in Eq. for particles with 
mv\/2 < U(r), z-motion is constrained within a V{t) 
puddle. It escapes when GC diffuses out of the U(r) 
puddle (Eig. Ic). Squeezing yields {vz)^ that vanishes 
and i;av ~ ('^gc)^ ^ Eq/B. As a result, Eq. ^ yields 
r< « ^A/{vgc)^ and 


DSC ^ 

XX 


2B ’ 


for Vz ^ (e 2 Vo/to)= v*, 


( 10 ) 


corresponding to electrons in Fig. lb with kz < kl; kl 
depends on the dispersion relation and v^. Eor 3DDM, 
hkl = {eTVQyjv^Y^. We note that Eq. can only be 
used for electrons with Vz ^ 'T’gc, Eq. (lb). However, in 
the opposite limit Vz fgc, ^z)\ is obviously smaller 
than ('i/gc)| in Eq. (j^, allowing us to neglect the former’s 
contribution, yielding as in Eq. (10). As a consis¬ 
tency check, we note that ^ for our regime of 
validity i > fc, eFo < ED- Therefore determines 
the range of electrons that obey Eq. (10). 

In contrast to Eq. (10), electrons with Vz ^ v^ do not 


have z-motion squeezed [case (b), see Eig. Id]. As a re¬ 
sult, GC samples many U(r) domains, with its x-y plane 
velocity scrambled over times r> ^ ^Jvz^ yielding an x-y 
plane mean free path I ^ v^ci/vz- This is captured in 
Eq. 0 whence (^c)c < giving Uav « {vz)(, ~ Vz- 

As a result, Eq. ^ yields oc IjB^. Importantly, for 
sufficiently large R, Vz > electrons, while mobile in 
the z direction, exhibit suppressed x-y plane mobility as 
compared with Vz <vl. As a result, class (a) trajectories 
dominates x-y plane transport. 

Linear Magnetoresistanee in 3DDM - To illustrate the 
striking effects of GC diffusion we specialize to 3DDM. 
Using the Einstein relation, and Eq. ( 8 pQ), we obtain 
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B 1 


( 11 ) 


where 5k is the electron energy, <a/(e^z/ 2 D) ~ 

Ro^- 4 /c*/( 47 r), and B « {Eq/vA) x hi{kF/kl). Here we 
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have used the 2D density of states for a kz slice in 3DDM 
as j/2d(m) = J2k^,ky '^(Sk - Ai) = M/27r?i^v|, and (A:^) 
obtained from the two trajectory classes (a) and (b). We 
note that the first term dominates over the second when 
B IB. It IS useful to re“VArite this condition as t'q 
where 1C = \/2(/i/eVb)^^^ x /k*])~^ ^ 1, since we 

are interested in /i > eVo- Hence, the first term always 
dominates in our regime of validity, ^ re- 

In the same way, the second term in Eq. 0 yields the 
usual expressions 

kz kz 


where ( = e^z/ 2 D(M)'^o/[ 2 (l + ^e have taken 

^cBr ^ 1 limit in the last expression. Since oc 

IjB^, for sufficiently large fields it provides a negligible 
contribution to axx as compared with Eq. 0 - 

An important diagnostic of magnetotransport is the 
Hall angle, tsuiOn = pxyjpxx = cfxylcfxx^ Using Eq. Q 
and writing a^x = (^fx (neglecting since uJcBv > 1 ), 
we obtain a 5-field magnitude independent 



where we have used n = /67r‘^h^Vp for a single fermion 

flavor in 3DDM, and gaussian correlated (U(r)U(r')). 
We note that the Hall angle changes sign when B field 
flips sign. Interestingly, the Hall angle can be tuned by 
Vq and p. Estimating Vb ~ 20 mV in 3DDM [ 22 ], and 
/i ^ 0.1 eV [T3IIT7] . we obtain tanOH ~ 2.4. 

We note that tunable Hall angle [Eq. 0 ] controls Q. 
Indeed, ^ is a non-monotonic function of Hall angle (and 
hence it depends on pjeVo), reaching a peak when Hall 
angle becomes unity. 

To summarize, we find that under the conditions 
^cBr ^ I 5 C ^ > eUo, G is independent of B 

field magnitude, yielding LMR according to Eq. 0 . Of 
the first two conditions, ^ ^ Tc can be expressed as 
B > Be = mvF/{e^) which is a more stringent condi¬ 
tion than uJcFtr ^ 1- This is seen by estimating Ttr using 
the Born approximation, giving Be that exceeds the field 
marking the onset of cJcTtr > 1 by a factor r\j {fi/eVo)\ 
Hence, we predict LMR as long as 5 > 5c. 

An important figure of merit for magnetoresistance is 
the ratio MR = [pxx{B) -Pxx{0)]/pxx{0)- Using Eq. 0, 
and noting the mobility 77 = cr 2 , 2 ,( 0 )/ne, we obtain 


MR = 


Pxx{B) pxx{^) ^ ( 77[cm^/Vs] \ ^ D[q^] X Q 

p,,(0) ^ V 104 ; ^ ^ 


- 20 X lO^cmVVs in 3DDM, Eq. ((M) 
5 — 100 at 5 = 10 T. Here we have 
used maximal Q = 1/2. Eor fixed 0, Eq. (14) is consistent 
with Kohler’s rule m since it scales with mobility. This 


Eor typical 77 ^ 1 
yields giant MR 


mirrors the experiments where MR scaled with temper¬ 
ature dependent mobility over wide range of 5 m- 

Another intriguing feature of LMR is an unconven¬ 
tional Hall resistivity at large fields. Conventionally, 
p^y = Bjne at large fields and is used to determine the 
density of carriers. However, since a^x oc 1/5 above, 
we obtain a Hall resistivity pxy = ^xyl{p\y + (^xx) ~ 

p^xyG 2 that differs from the conventional case by a fac¬ 
tor G 2 = tan^H X G- Since ^2 < 1, small Hall angles, 
can lead to unconventional pxy. Indeed, in these situ¬ 
ations, determining density through Hall measurements 
via n = B/{pxye) may overestimate the density. 

We note that for other dispersions, our treatment 
above follows through yielding LMR under the same con¬ 
ditions of uOeTtr ^ 1, ^ Tc, /i > eVb- Howcvcr, features 

e.g., /c*, a and Hall angle are altered appropriately. 

Inelastic scattering - Energy relaxation through in¬ 
elastic scattering (e.g. through phonon scattering) in the 
2 )-direction can drastically affect GC motion by mixing 
squeezed ^ with unconstrained 2 : trajectories. Absorp¬ 
tion of phonons with Tiuj — ^v‘1 relaxes the energy 
constraint [Eq. 0 ], allowing Vz < electrons to jump 
out of U(r) troughs. If r< « ^/ugc > ri„(e = Vq - 
these electrons exhibit Dxx{kz) oc Here ri„(e) is 

the time for an electron to absorb energy 5 . Suppressed 
at low T, phonon-assisted escape leads to LMR degrada¬ 
tion when ksT ^ Vb- 

However, when typical Vb ^ phonon-assisted es¬ 
cape becomes difficult even at high temperatures since 
the maximum energy that can be absorbed from phonons 
is fico’D ^ Vb; cjd is the Debye frequency. As a result, in 
this regime LMR is stable even at high temperatures and 
large fields, as recently observed in Cd 3 As 2 mills] where 
LMR persisted at 300 K and high fields. 

While detrimental to LMR from GC diffusion in 3D, in¬ 
elastic scattering can have the opposite effect in 2D. Con¬ 
ventionally in 2D, GCs form closed orbits along equipo- 
tential lines of a disorder potential yielding localization 
behavior [28]. However, for inelastic phonon scattering 
which is not so strong to entirely disrupt the GC mo¬ 
tion, but strong enough {^/vge ^ ^in) to induce switch¬ 
ing between adjacent equipotentials [29], the GC trajec¬ 
tories can become open, moving through multiple U(r) 
domains. In this regime, we speculate 5^^ ^ Vge^ oc 1/5 
as above. Interestingly, 2D semi-classical regime LMR 
was reported previously [30] . 

Semi-classical GC diffusion can conspire to produce 
LMR in metals. Importantly, the requirements for GC 
magnetoresistance are modest - arising in the semi- 
classical regime with multiple occupied Landau levels, 
and for weak and smooth disorder. Giant MR ratios, 
5-field magnitude independent Hall angles, p and Vb 
tunability, and stability at high temperatures, make 
GC diffusion and its magnetoresistance easy to iden¬ 
tify in experiment. Indeed, these features bear strik- 
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ing resemblance to LMR measured recently in a vari¬ 
ety of 3DDM [TSHTT]. Additionally, oscillatory motion 
of the GC trajectories along B could have interesting, 
polarization-dependent, absorption signatures in the Ter¬ 
ahertz regime. 
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